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Solution to the Volterra Matrix Equation 
of the 1st kind with Piecewise Continuous Kernels 

Denis Sidorov 



In this paper we address the following matrix Volterra integral equation 
^ ■ of the first kind 



t 



J K(t, s)x(s)ds = f(t), < t < T. (O.i; 



c/r o 

We assume mxm matrix kernel K(t,s) has jump discontinuities on curves 



s = oti(t), i = 1, n — 1, < s < t < T as follows 



in- K{t,s)=< 
m 

o 

(N 



Ki(t,s), 0<s<ai(t), 
K 2 (t, s), ai(t) < s < a 2 {t), 



(0.2) 



K n (t, s), a n _i(t) < s < t, 



/(*) = ■ • ■ , /mW)' a?W = (a^i W, • • • , Matrices ^(t, s) are 

mxm. We assume that these matrices are defined, continuous and have 
continuous derivatives wrt t correspondingly in Di = |s, t\cti-i(t) < s < 
ai(t)\, i = l,n, ao = 0, a n (t) = t. Functions fi(t),ai(t) have continuous 
derivatives, / f (0) = 0, ^(0) = 0, < a[{0) < a 2 (0) < ■■■ < <_i(0) < 1, 
< OLi{t) < a 2 (t) < ■■• < On-lit) < t for t G (0, T}. In this chapter 
we search for the continuous solutions of the equation ( 10.11 ) for t G (0, T'], 
where < T' < T, lim x(t) can also be infinite. We assume below that 
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1001-012. 
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each matrix Ki(t, s), i = 1, n has continuously differentiable wrt t exten- 
sion on the interval < s < t < T. Homogeneous systems can have 
nontrivial solutions. Differentiation of the system ( 1DTTT) wrt t in contrast 
to the classic case 0,E],[E] lead us to the new class of the Volterra integral 
equations with functionally perturbed argument. That is the reason why 
the conventional techniques are not applicable for the system ( 10.11 ) and 
studies of the systems with jump discontinuous kernels ( 10.21 ) are theoreti- 
cally and practically interesting object to study. Operator equations with 
functionally perturbed argument are addressed in [TTl[T2"]. 

The paper is organized as follows. In section 1 we propose the algorithm 
for construction of the logarithmic power asymptotics 

N 

x(t) = ^^(lnt)f (0.3) 

i=0 

of the desired continuous solutions of the system ( 10.11 ) . In the section 2 we 
proof theorem of existance of the parametric families of solutions to the 
system ( 10.11 ). Finally we derive the sufficient conditions for existance and 
uniqueness of continuous solution. 

1 Construction of asymptotic approximation 
of solutions to inhomogeneous system 

Let the following condition be fuilfilled 

N 

(A) Matrices V%= Y Ki Vfl t v , i = 1, n exist, vector-function / (t) = 

N N 

f v t v , and polynomials af(t) = Yl a t^"i i = 1, n — 1, where < 
an < Q2i < • • • < a n -i,i < 1, such as for t — > +0, s — > +0 the 
following estimates are fulfilled Ki(t, s)—Vi(t, s) = 0((t J rs) N+1 ), i = 
T^, f(t) - f N (t) = (D(t N+1 ), ai(t) - af{t) = (D(t N+1 ), i = T~n^l. 
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Expansions in powers of t, s in (A) we call below as "Taylor polynomi 
als" of the corresponding elements. Let us introduce the matrix 

ra-l 

B(j) = K n (0 : 0) + J2(®'M) 1+J ( K ^ 0) - K i+1 (0 : 0)) (1.1 

i=l 

and the following algebraic equation 

L(j) d ^detB(j) = 0. 

We call it as characteristic equation of the system of integral equations 
( 10.1 1 ). Since /(0) = matrices Ki(t,s) and vector f(t) have continuous 
derivatives wrt t, then differentiation of the system ( 10.11) lead us to the 
equivalent system of integral-functional equations 

71—1 

F(x) = K n (t,t)x(t) + Y,a' i (t){K i (t,a i (t))- (1.2) 

i=l 
n a f 

-K i+1 (t,ai(t))}x(ai(t)) + Yl J K$ 1 \t 1 s)x(s)ds-f(t) = 0, 

where cvq = 0, a n (t) = t. 

Here we do not assume that homogenius system of ( 10.11 ). has only 
the trivial solution. Hence the homogeneous integral-functional system 
corresponding to ( 11.21) can have nontrivial solutions. Here we follow [TTl 



and seek an asymptotic approximation of a particular solution of the 
inhomogeneous equation ( 11.21) as following polynomial 

N 

x (t) = j2 x ji int ) tJ - ( L3 ) 

Let us demonstrate that coefficients Xj depend on In t and free parameters 
in general irregular case. That is in line with possibility of nontrivial 
solutions existance of the homogeneous system. 



The regular and irregular cases are possible when compute the coeffi- 
cients Xj. 

1.1 Regular case: L(j) ^ 0, j e (0, 1, . . . , N) 

In this case the coefficients Xj are constant vectors from M. m . Indeed, let us 
substitute an expansion ( 11.31 ) in the system ( 11.21 ) and apply the method of 
undetermined coefficients taking into account the condition (A). Then we 
get the recurrent sequence of linear system of algebraic equations (SLAE) 
with respect to the vectors Xj : 

B(0)x = f(0) J (1.4) 

B(j) Xj = Mj(x 0: . . . , Xj -!), j = 1, . . . , N. (1.5) 

Vector Mj can be expressed in a certain way through the solutions 
Xq : ... : Xj-i of the previous systems and coefficients of the "Taylor poly- 
nomials" from the condition (A). 

Since in the regular case det B(j) ^ 0, then vectors Xq : . . . , xn can be 
uniquely determined and asymptotic ( 11.31 ) can be therefore constructed. 

1.2 Irregular case: equation L(j) = has integer 
roots. 

Let us introduce the definitions: 

Definition 1. We call quantity j* the regular point of the matrix B(j) if 
matrix B(j*) is invertible. 

Definition 2. We call quantity j* simple singular point of the matrix 
B{j) if det SOT) = °> det[(£%-*)^)]^ =1 ^ 0, where {^}\ is basis 
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in N(B'(j*)), B'(j*) is transposed matrix, and B^\j) is derivative of the 
matrix wrt j. 

Definition 3. We call quantity j* k + 1-multiple singular point of the 
matrix B(j) if det B(j*) = 0, derivatives B^ l > (j *),..., B^ k \j*) are zero 
matrices, 



det 



(£ (fc+1) Cf)<^« 



r 



i,Z=l 



/c > 1, is basis in N(B(j*)), {^K is basis in N(B'(j*)). 

n-1 

Let us notice that 5 (fc) (j) = E °) ~ ^+i(0,0)), 

where = lnc^(0). 

Remark 1.1. Let m = 1 (this case correspond to the single equation ( 10.11 )). 
Then = Hence in case of the single equation the definition 2 

means that j is the single root of the characteristic equation L(j) = 0, 
and definition 3 means that j is (k + 1) multimple root of this equation. 

Let us demonstrate that in irregular case the coefficients Xj are poly- 
nomials in power of Int and depend on arbitrary constants. Order of poly- 
nomials and number of arbitrary constants are associated with muliplicity 
of singular points of the matrices B(j) and ranks of these matrices. 

Indeed, since the coefficient xq in irregular case may depend on lnt, 
then we may employ the method of undetermined coefficients and seek xq 
as solution of the difference system 

n-1 

tf n (0, 0)x (z) + ^ «K0)(^(0, 0) - K i+1 (0, 0))x (z + <*) = /'(0), (1.6) 

i=i 

where di = lna'(0),z = Int. 

Here there are three cases: 
1st Case. 

L(0) 7^ 0, i.e. det 5(0) ^ 0. Then the coefficient xq does not depend 
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onz and determined uniquely from the SLAE ( 11.41) with invertable matrix 

B(0). 

2nd Case. 

Let j = be the simple singular point of the matrix B(j). We seek the 
coefficient xq(z) from the system ( 11.61) as linear vector- function 

x (z) = x 01 z + x 02 . (1.7) 



Let us substitute ( 11.71) in ( 11.61) . and we get the following two SLAE for 



determination of the vectors xqi, X02'. 

B(0)x 01 = 0, (1.8) 

B(0)x 02 + B^(0)x 01 = f'(0). (1.9) 

r 

Here detB(O) = 0, {4>i} r i is basis in N(B(0)). Hence xo\ = ^ c fc0fc- 

k=i 

Vector c = (ci, . . . , c r )' can be determined uniquely from the conditions of 



the system ( 11.91) resolvability, i.e. fro the SLAE 

r 

^(B«(O)0 fe , ^)c k = (/'(0), * = 
fe=i 

with non-singular matrix. The coefficient X02 can be determined from the 
system ( 11.91) with accuracy up to span(0i, . . . , <fi r ). Therefore, in 2nd case 
the coefficient xo(z) will depend on r arbitrary constants and linear wrt z. 

3rd Case. Let j = be the singular point of the matrix B(j) of multi- 
plicity k + 1, where k > 1. Solution £q(z) to the difference system ( 11.61) 
we search as polynomial 

X (z) = X iZ k+1 + X 02 Z k H h X 0k+ iZ + (1-10) 



6 



Let us substitute the polynomial ( 11.101 ) into the system ( 11.61 ). We next 
take into account the equality 

— k B(j) = ^(a;.(0)) 1+ X fc (^(0,0) - K m (0,0)), 



dj 



i=l 



where = In a^(0) and equation coefficients of powers z +1 , z fc , . . . , z, z° to 
zero yields the recursive sequence of linear algebraic equation with respect 
to the coefficients #01, ^02, • • • , ^ofc+2 : 
B(0)x Q1 = 0, 

k + 1 



5(0)x 02 + 5( 1 )(0)(^ 
B(0)a;oi+i + 5(0(0) 
... + 5«(0) 



fc + 1 
fc + l-Z 
fc+1 -Z+l 
+ 1 -Z 



x i = 0, 



x 01 



+ 5^(0) 



fc 

fc + l-Z 



£02 + 



x i — 0, / — 



B(0)x 0fc+2 + 5( fc+1 )(0)x l + 5^X02 + • • • 5W(0)xo fc+ i = /'(0). (1.11 



d l B(j) 
df 



Since in this case due to the conditions of the definition 3 the derivatives 
, i = are zero matrices, then 



j=0 



x j = ^ Cfj^j-, i = 1, . . . , + 1. 
i=i 



As result the system ( 11.111 ) can be rewritten as 



1.12) 



Since det 



;5( fc+1 )(o)0^ fc ) 



7^ 0, then vector c 1 = f (en, . . . , ci r )' can 



i,k=l,r 



be uniquly determined from the conditions of resolvability of the system 
021) . Therefore, 

r 

X0k+2 = ^ Cfc+2 j<ftj + X k+ 2, 
.7=1 
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Xk+2 is the particular solution to the SLAE ( 11.121 ). Vector c k+2 = f 
(cfc+2,1, • • • Ck+2,r)'i an( i vectors c 2 , . . . , c k+1 , remain arbitrary. Therefore 
in 3rd case the coefficient Xq(z) is polynomial of power of k + 1 wrt z and 
it depends on r x (h + 1) arbitrary constants. One may apply the method 
of undetermined coefficients and take into account the following equality 



/k 
t j \n k tdt = t j+1 ^2(-l 
s=0 



,k(k-l).-.(k-(s 



u + 1: 



s+l 



i»ln*- t . 



and construct the system of difference equations for determination of the 
coefficent x\{z) [z = hit). Indeed, 



Fix) 



dcf 



x=Xo+x%t 



n-1 



ir n (0,0)x!(z) + ^(a;(0)) 2 (^(0,0)- (1.13) 



i=i 



K i+1 {0, Q))x!(z + a,) + Pi(x (z)) 



t + r(t), r{t) =o(t). 



Here P\(xq(z)) is certain polynomial on z, which power is equal to the 
multiplicity of the singular point j = of the matrix B(j). From ( 11.131 ) 
due to the estimate r(t) = o(t) for t — > it follows that the coefficient 
x\(z) must satisfy the the following system of difference equations 



n-1 

tf„(0, 0)^(2) + ^(a'(0)) 2 (^(0,0)-^ +1 (0,0))x 1 (z + a l )+ (1-14) 

+Pi(x (z)) =0. 

If j = 1 is regular point of the matrix B(J), then system ( 11.141 ) has solution 
x\(z) as polynomial of the same order as multiplicity of the singular point 
j = of the matrix B(0). If j = 1 is singular point of the matrix B(j), 
then solution xi(z) can be constructed as polynomial of power ho + hi, 
where ho and k\ are multiplicity of singular point j ' = and j = 1 of the 
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matrix B(j) correspondingly. The coefficient x\(z) depends on ro&o + ^i^i 
arbitrary constants, where tq = dimiV(.B(0)), T\ = dimN(B(l)). 
Let us impose the following condition 

(B) Let matrix B(j) has only the regular points in the array (0, 1, ... , N) 
or has singular points of multiplicities kj. 

Then, in a similar way we can calculate the remaining coefficients 
3^2(^)5 • • • 5 %n{ z ) of x(t) from the sequence of difference equations 

71—1 

K n (0, 0) Xj (z) + ^(c/(0)) 1+J (^(0, 0) - K i+1 (0, 0)) Xj (z + a,)+ 
i=i 

+Vj(x (z), . . . , ^-1(2;))) = 0, j = XN. 
Hence we have the following 

Lemma 1.2. Let conditions (A) ; (B) be fuilfilled. Then exists vector 

N 

function x(t) = ^2 Xi(lnt)t\ such as for t — > +0 the following estimate is 

satisfied \F(x(t))\^m = o(t N ) and the coefficients Xi(lnt) are polynomials 

of Int with insreasing powers smaller then sum kj of the singular points 

j 

j of the matrix B(j) from the array (0, 1, . . . , i). Coefficients Xi(lnt) de- 

i 

pend on ^ dimN(B(j))kj arbitrary constants. 
3=0 

2 Existence of continuous parameter family 
of solutions theorem 

Since < a-(0) < 1, ^(0) = 0, i = 1, n — 1, then for any < e < 1 exists 
T' £ (0, T] such as ma x I^COI — £ an< ^ SU P ^t^~ — £ - 

i=l,n-l,te[0,T'} i=T^n-L,te(0,T'] 

Let the following condition be fulfilled 
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(C) Let det K n (t,t) ^ 0, t G [0,T / ] and N* is selected large enough to 
have the following inequality be satisfied 

max £ JV *|i^~ 1 (t,t)|£(M TO ^ Mm )X (2.1) 
te[o,T] 

n-l 
i=l 

where | • \c(u.m^ m ^ is norm of the m x m matrix. 

Lemma 2.1. Le£ t/ie condition (C) 6e fulfilled. Let C(o,T'l fre c/ass of 
the vector functions continuous for t G (0, T"] and and having a limit 
(possibly infinite) int — > +0. Let in C(o,t'1 exists the element x(t) such as 
the following estimate 

\F{x{t))\ Rm = {t N ),N>N* 
is fulfilled. Then equation ( 11.21) in C(q,t'] has the solution 



x{t) = x{t) + t N *u{t), (2.2) 

where u(t) G Cr 0) T'i and can be uniquely defined with successive approxi- 
mations. 

Proof. Let us substitute ( 12.21 ) in the equation ( 11.21) . We get the following 
integral-functional system for determination of u{t) 

K n (t,t)u(t) + ^a[(t)(^Y (k^o*®)- (2.3) 



i=i 



■Ki+i(t,ai(t)) )u(ai(t))+ 



+ J2 J K ^^ s \y <s)ds + F(x{t))/t N * = 0. 



ai-i(t) 
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Let us introduce the linear opearators: 

n-l / s.w N* 



Lu A M K -% t) £ a[{t) {^Pl <*(*))- 
-Ki+i(t,ai(t))ju(ai(t)), 

Oli(t) 

,_i ai_i(t) 

Then we can present the system ( 12.31 ) as follows 

u + (L + K)u = 7(f), 

where 7(t) = K r ^ 1 (t,t)F(x N \t))/t is continuous vector function. Let us 
introduce banach space X of continuous vector-functions u(t) with norm 

\\u\\i = max e~ lt \u(t)\w™, I > 0. 

0<i<T' 

Hence because of inequalities sup < £ < 1 and condition C) for 

te[o,T] 

V7 > norm of linear functional operator L satisfies the following estimate 

\\L\\c(x^x) <q<l- 

Moreover for integral operator K for large enough / the following estimate 
is fulfilled 

\\K\\c{x^x) <qi<i-q- 

As result for large enough / > we have 

\\L + K\\ L{X ^ X ) < 1, 

i.e. the linear operator L + K is contracting in X. Hence the sequence 
u n = — (L + K)u n -\ + j(t), uq = j(t) converges. 

□ 
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Theorem 2.2. (Main Theorem). Let the following conditions (A) ; (B) ; 
(C) ; /(0) = be fulfilled. Let also the matrix B(j) has exactly v singular 



points ji, . . . , j v of multiplicities ki, i = in the array (0, 1, ... , N)) 
and the rest of the points of the array are regular. Let rankB(ji) = 



i = 1, v. 

Then equation ( 10.11 ) for < t < T' < T has he solution 



x (t) = x(t) +t N *u(t), 

V 

depending on ^(m — rf)ki arbitrary constants. 

Proof. Due to the conditions of the Lemma 1 and imposed conditions of the 
theorem it make it possible the construction of asymptotic approximation 
x(t) of the desired solution in the from of logarithmic power polynomial 

N 

Xi(lnt)t\ Moreover the coefficients Xj(lnt) will depend on the specified 

i=0 

number of arbitrary constants. Because of the Lemma 2 we can apply the 
substitution x(t) = x(t) +t N *u(t), and continuous function u(t) can be 
constructed with successive approximations. 

□ 

Remark 2.3. In the main theorem's conditions for the asymptotic approx- 
imation x{t) of the desired solution the following asymptotic estimate 
\x{t) - x{t)\ Rn = 0{t N *), +0 is fulfilled. 



Corollary 2.4. Let ai(t) = a^t, i = 1, n — 1, < a\ < «2 < ••• < 



a n -i < 1, elements of the matrices Ki(t, s), i = 1, n — 1 and vector func- 
tion f(t) has an analytic extension to \s\ < T, \t\ < T, /(0) = 0. Matrix 
K n {t,ty l is analytical for \t\ < T. Let det B(j) ^Oj'eNUO. Then 



equation ( 10.11 ) has unique solution x(t) = ^ xft for < t < T. 

i=0 

In some cases, the conditions of Theorem 1 we can construct a para- 
metric family of solutions in closed form. 
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t/2 t 

Example. Let us consider the system J Kx(s)ds + J (K — 2E)x(s)ds = 

t/2 

dt, < t < oo, where K is symmetric constant matrix m x m, d 6 
M. m , x{t) = (xi(t), x m (t))', 1 is eigenvalue of the matrix K of rank 
r, {0i,...,0 r } is corresponding orthonormal system of eigenvectors. This 
system has parametric family of solutions 



> (d c 

ar(t) = - hit } / ' t + ci0i + ... + c r r + d. 
L — ' In z 

i=i 



Here C\ . ... , are arbitrary constants, vector a satisfies the SLAE 

r 

(K - E)a = d-^2(d : 



i=l 



3 Sufficient conditions for existence of a 
unique continuous solution 

For sake of clarity in this section let us suppose ai{t) = a{t, i = 



1, n — 1, < a,\ < oii < ■■■ < OL n -\ < 1. Let us introduce the matrix 

71—1 

D(t) = ^2a i K-\t : t)(K i (t : a 1 t) - K l+1 {t, ai t)). 

i=l 

Let the following condition be fulfilled: 

(S) \D(0)\ C (Rm^ R m) < 1; sup \K-\t,t)K(t,s)\£(Km-Mm) < c < oo. 

0<s<t<T 



Here and below the matrix K(t, s) is defined with formula ( 10.21) 



Theorem 3.1. (Sufficient conditions for the existence and uniqueness of 
solutions). Let the conditions (S) be fulfilled, all the matrices Ki(t : s) are 
continuous in ( 10.21) and they are continuously differentiable wrt t, vector 
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f(t) has continuous derivative, f(0) = 0. Then an equation ( 10.11) has 
unique solution in CrnTl- Moreover, the solution can be found by the method 
of steps, combining it with the method of successive approximations. 



Proof. Let us rewrite the equation ( 11.21) which equivalent to the equation 
( 10.11) as follows 

x{t) + Ax + Kx = J{t), (3.1) 
where the following notations are introduced 

def 

Ax = K~ l {t,t) ai(Ki(t,ait) - Ki+i(t,ait))x(ait), 
i=i 

^ = E / K-%t)K^\t : s)x(s)ds : f(t) = K-\t,t)f^(t). 

i=lcti—it 

Let us fix q < 1 and select hi > such as max |-D(^)l£( ffim ->- Rm ) = 
q < 1. Because of the condition (S) exists such h\ > 0. Let < h < 
min{/i!, ^^}, where constant c is defined in the condition (S). Lets divide 
the interval [0, T] into the intervals 

[0,h], [h,h + eh], [h + eh, h + 2eh] , . . . (3.2) 

where e is selected from (0, 1] such as a n -\ < jj-. Denote by xo(t) the 
restriction of the desired solution x(t) on the interval [0, h] and let x n {t) 
be the restrictions to intervals 

I n = [(1 + (n - l)e)h, (1 + ne)h] , n = 1, 2, . . . . 

Because of the selection of e for t £ I n the "perturbed" argument a{t G 

71-1 

(J Ik- Such an inclusion makes it possible to apply the known method of 

k=i 
steps 



For construction of the element xo(t) G Crow we construct the sequence 
KM) = 

x n (t) = -Ax%- 1 -Kx n - 1 + J(t), 
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x° Q (t) = f(t),te[o,h]. 

Because of the selection of h we have an estimate 1 | \c(c [Q h] ^c [0 h] ) < 1- 

Hence exists the unique solution Xa(t) of the equation ( 13.1 1 ) for t G 

[0,/i]. Sequence Xq^) converges uniformly to that unique solution xo(t). 

Let us continue the process for the desired solution construction for t>h, 

i.e. on the intervals I n , n = 1, 2, . . . . To be specific let e = 1 in ( 13.21) . 
Once we get a?o(£) £ Qo w computed we constract a?i(t) in the space 

^[/i,2/i]- We find from the Volterra integral equation of the 2nd kind 
t h 

x(t)+ J K~ 1 (t J t)K' t (t,s)x(s)ds = J(t)-Ax - J K- l {t,t)K[{t, s)x (s)ds 
h o 

using successive approximations. In this case xo(h) = x\{h). 
Let us introduce the continuous function 

|,„W, o<t<k, 

{ Xl (t), h<t<2h, 
which is Is the restriction of the desired continuous solution x(t) on the 
interval [0,2/i]. Then element X2(t) G CW^i can be computed with suc- 
cessive approximations from the Volterra integral equation of the second 
kind 

t 2h 

K- 1 (t,t)K' t (t,s)x 1 (s)as 

2h 

Finally we constract the desired solution x(t) G C^t] to the equation 

Z' 
h 



i zn 

x(t)+ J K r ^ 1 (t,t)K' t (t,s)x(s)ds = J(t)-Ax 1 - J K~\t,t)K' t (t, s) Xl (s)di 



( ICUl ) during the N steps (N > f ). 



□ 



Example. Integral equation 

t/2 t 



J Kxit - s)x{s)ds + J K 2 {t - s)x(s)ds = f(t), < t < T, 



t/2 
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where K\{t — s) = K^it — s) + E, are matrices m x m, E is unit matrix, 
l-^-2~ (0)1 £(R m ->-K m ) < 2, matrix K2(t) and vector function /(£) have con- 
tinuous derivatives wrt t, f(0) = 0, satisfies the Theorem 2 conditions and 
has the unique continuous solution. 
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